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A  COMPARISON  OF  TWO  TECHNIQUES  FOR  THE 
INTERPOLATION/EXTRAPOLATION  OF  FREQUENCY  DOMAIN  RESPONSES 


Sharath  Narayana,  Tapan  K.  Sarkar  and  Raviraj  Adve 
Department  of  Electrical  and  Computer  Engineering 
121  Link  Hall 
Syracuse  University 
Syracuse,  New  York  13244-1240 


ABSTRACT:  In  this  paper,  a  comparison  is  made  between  two  methods  which  are  used 

for  interpolation/extrapolation  of  frequency  domain  responses. 

The  first  is  the  direct  method  based  on  the  principle  of  a  model  based  parameter 
estimation  and  the  second  method  uses  the  properties  of  discrete  Hilbert  transfomis  based  on  the 
principle  of  causality. 

1.  INTRODUCTION 

In  a  host  of  applications  in  engineering,  it  is  necessary  to  obtain  infoimation  about  a 
system  over  a  broad  range.  In  most  cases,  it  is  not  possible  to  evaluate  the  parameter  of  interest 
in  a  closed  form.  However,  either  theoretical  or  experimental  data  is  available  in  a  narrow  band. 
The  first  method  uses  the  principle  of  analytic  continuation  to  extrapolate/interpolate  the  data 
over  a  wide  band.  The  second  method  uses  the  property  of  Hilbert  transfonn  which  relates  the 
real  and  imaginary  components  of  the  frequency  domain. 

In  Section  2,  the  direct  method  is  studied.  Section  3  describes  the  numerical 
implementation  of  the  Hilbert  transform.  Section  4  describes  the  algorithm  used  to  implement 
this.  In  Section  5,  few  numerical  examples  are  studied  and  comparisons  made. 
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2. 


THE  CAUCHY  METHOD 

Consider  a  system  function  H(s).  The  objective  is  to  approximate  H(s)  by  a  ratio  of  two 
polynomials  A(s)  and  B(s)  so  that  H(s)  can  be  represented  by  fewer  variables. 

Hence,  consider 


H(s) 


A(s)  _  ^=0 


B(s) 


b  s 

■^k=0 


(1) 


Here,  the  given  information  could  be  the  value  of  the  function  H(s)  and  its  Nj  derivatives  at  some 
frequency  points  Sj,  j=l,  ...J.  If  H"(Sj)  represents  the  n”*  derivative  of  H(s)  at  point  s=Sj,  the 
Cauchy  problem  is: 

Given  H<'’>(Sj)  for  n=0,  ...Nj,  J=l,  ...J,  find  P,  Q,  {a„  k=0,  ...P},  and  {b^,  k=(),  ...Q). 

The  solution  for  {aj.}  and  (bj  is  unique  if  the  total  number  of  samples  is  greater  than  or 
equal  to  the  total  number  of  unknown  coefficients  P  -  Q  +  2  [Ij,  i.e. 

J 

N  =  E  (R  + 1)  i  P  +  Q  +  2 

J=i 

By  enforcing  the  equality  in  equation  (1)  one  obtains 

A(s)  =  H(s)B(s) 

Differentiating  the  above  equation  n  times,  and  evaluating  the  expressions  at  point  Sj, 
results  in  the  binomial  expansion, 

n 

A(“^(s.)  =  E°CH("-‘>(s.)B®(Sj) 

^  i=0 

where, 


n! 


n!  represents  the  factorial  of  n. 

Using  the  polynomial  expansions  for  A(s)  and  B(s),  equation  (3)  can  be  written  as 


2 


(4) 


where, 


p 

E 

k=0 


‘■(j.n),k*k 


Q 

=  E 

k=0 


®(j.n),k'>k 


A  = 

Sj^  u(k-n) 

(5) 

(k-n)! 

A.«.v  =  ^  C,  (Sj)  u(k-i)  ^  s<‘-l> 

i=o  ^  (k-i) ! 

(6) 

n=0,  1,  ...Nj,  j=l,  ...J,  where  u(k)=0  for  k<0  and  = 

Define, 

=  1  otherwise. 

[■^a.n).0’ 

(7) 

(8) 

[a]  =  [ao,aj. 

a2,...apf 

(9) 

[b]  =  [bo,bj,b2,...bQf 

(10) 

The  order  of  matrix  A  is  N  x  (P  +  1)  and  that  of  B  is  N  x  (Q  +  1).  Then,  equation 

(4)  becomes 

[A]a  =  IB]b 

(11) 

or 

[Al-Bl  j 

1 

=  0 

) 

(12) 

For  ease  of  notation,  define  [C]  =  [A]  -  B],  C  is  of  order  N  x  (P  +  Q  +  2).  A  Singular  Value 
Decomposition  (SVD)  of  the  matrix  C  will  give  us  a  gauge  of  the  required  values  of  P  and  Q 

[2],  A  SVD  results  in  the  equation 
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IUI[E1[V]« 


a 


(13) 


b 


=  0 


The  matrices  U  and  V  are  unitary  matrices  and  Z  is  a  diagonal  matrix  with  the  singular 
values  of  C  in  descending  order  as  its  entries.  The  columns  of  U  are  the  lelt  eigenvectors  of  B 
or  the  eigenvectors  of  CC".  The  columns  of  U  are  the  right  eigenvectors  of  C  or  the 
eigenvectors  of  C"C.  The  singular  values  are  the  square  roots  of  the  eigenvalues  of  the  matrix 
C"C.  Therefore,  the  singular  values  of  any  matrix  are  real  and  positive.  The  number  of  nonzero 
singular  values  is  the  rank  of  the  matrix  in  equation  (12)  and  so  gives  us  an  idea  of  the 
information  in  this  system  of  simultaneous  equations.  If  R  is  the  number  of  nonzero  singular 
values,  the  dimension  of  the  right  null  space  of  C  is  P  +  Q  +  2  -  R.  Our  solution  vector  belongs 
to  this  null  space.  Hence  to  make  this  solution  unique,  we  need  to  make  the  dimension  of  this 
null  space  1  so  that  only  one  vector  defines  this  space.  Hence  P  and  Q  must  satisfy  the  relation 

R  +  l=  P+  Q+  2  (14) 


The  solution  algorithm  must  include  a  method  to  estimate  R.  This  is  done  by  starting  out 
with  the  choices  of  P  and  Q  that  are  higher  in  dimension  than  can  be  expected  for  the  system  at 
hand.  Then  we  get  an  estimate  for  R  from  the  number  of  non-zero  singuhir  vtilues  of  the  matrix 
C.  Now,  using  equation  (14)  we  get  better  estimates  for  P  and  Q.  Letting  P  and  Q  stand  for 
these  new  estimates  of  the  polynomial  orders,  we  can  recalculate  the  matrices  A  and  B. 
Therefore,  we  come  back  to  the  relation 


MA|  -B] 


=  0 


(15) 


IC]  is  a  rectangular  matrix  with  more  rows  than  columns.  Many  methods  to  solve  equation  (15) 
are  well  documented  [2].  For  reasons  indicated  in  the  appendix,  we  choose  the  method  of  Total 
Least  Squares  (TLS)  [3]. 


3  TRANSFORM  RELATIONSHIPS 

This  section  briefly  covers  some  of  the  properties  of  sequences  and  their  Fourier 
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transforms.  Any  complex  sequence  h[n|  can  be  expressed  as  a  sum  of  a  symmetric  sequence 
hjn]  and  an  anti-symmetric  sequence  h„|  n].  In  the  case  of  real  sequences,  these  are  called  even 
and  odd  sequences  [4],  Therefore, 

h[n]  =  h^[nl  +  h„[n]  (16) 

hjn]  =  hj-n]  (17) 

hjn]  =  -hj-n]  (18) 

The  Fourier  transform  of  any  complex  sequence  hfn]  is  represented  by  H(e'“),  where 

00 

H(el“)  =  X)  h[n]  6^“"  (19) 

n=-oo 


Therefore, 


This  implies, 


and 


00 

H(e^“)  =  h[n] 
n=-<» 


H«(ej“)  = 


(20) 


(21) 


jHg(ei“)  =  (22) 

Also,  for  real  h[n],  Hg,(e’“)  =  H3i(e  j“)  which  is  an  even  function,  and  H3(e'“)  =  -113(0^“)  which 
is  an  odd  function. 


3.1.  THE  HILBERT  TRANSFORM  RELATIONSHIP  FOR  A  CAUSAL  SEQUENCE 
Consider  a  periodic,  real,  time  domain  sequence  hp[n]  with  period  N,  that  is  related  to  a 
finite  length  sequence  h[nl  of  length  N  by 
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00 


(23) 


hp[n]  =  h[n+iisri 

i=-oo 

From  Eq.  (16)  we  have, 

hp[n]  =  hp,[n]  +  hp^[n] 

and  from  Eq.  (17)  we  have 

VW  =  ^(hp[n]  +  hp[-n]} 

and  similarly  from  Eq.  (18), 

lyM  =  -  V-xll 

If  we  have  N=2r  (where  r  is  a  positive  integer)  then, 

2hp,[nl  n  =  1,  ...  |-1 
hp[n]  =  hpeM  "  0’  f 

0  n  =  -  +  1,  ...  N-1 

2 

Also  the  odd  part  of  the  sequence  can  be  expressed  as, 

h^[n]  n  .  1,  ...  |-1 

„  n  N 

vw  H  0  “  =  0-  y 

-ly[nl  n  =  ^  +  1.  ...  N-1 

If  we  define  UpN[n]  as  the  periodic  sequence 


(24) 


(25) 


(26) 


(27) 


(28) 
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UnN[n]  =  i  1 


n  =  1,  — -1 

2 

n  N 

n  =  0,  — 

2 


n  =  -  +  1,  ...N-1 
2 


we  can  express  h  [n]  as 


hp[n]  =  hpeMUp^Ln] 


Equivalently,  the  Fourier  transform  of  h  [n]  yields, 


IN-i 

“  TT  “  ^p«M  +  jHpg[k] 

m=0 

Eq.  (29)  can  be  alternatively  expressed  as, 

u^[n]  =  2u[n]  -  2u[n-^]  -  6[n]  +  6[n--l 

2  2 

where  u[n]  is  the  unit  step  sequence,  and  5[n]  is  the  unit  sample  sequence  1 10]. 
The  Fourier  transform  of  UpN[n]  is. 


-2jcot  —  k  odd 

=  1  N  J 


k  even 


with  UpN[0]  =  N  as  derived  from  the  definition  of  UpN[k]. 
Defining 


-  N6(k) 


we  have,  from  Eq.  (31) 
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Hp[k] 


N-1 

Hp«[k]  ^  ^  E 

^  01=0^5^]^ 


(35) 


Hence  we  get  the  imaginai^  part  of  the  Fourier  transform  as, 

N-1 

jHpsiM  =  i  E  (36) 

N  m=0,m^k 

This  is  the  Hilbert  transform  relationship  between  the  real  and  imaginary  parts  of  the 
Fourier  transform  of  a  periodically  causal  sequence.  If  h[n]  =  0  for  n  <  0  and  for  n  >  N/2  then 
the  periodicity  maybe  removed  and  we  have  the  same  relationships  between  the  real  and 
imaginary  parts  of  the  Fourier  transform  of  h[n]. 


.  H-[m]VN[k-m]  0  ^  k  ^  N  -  1 

jH  [k]  =  Jn  J  NL  (37) 

[  0  otherwise 

In  the  previous  expression  the  cotangent  term  is  itself  periodic  with  period  N,  so  when 
we  compute  the  N  point  discrete  Fourier  transform  of  the  real  causal  time  sequence,  the 
relationship  between  the  real  and  imaginary  parts  will  be  affected.  Eq.  (37)  is  the  Hilbert 
transform  relationship  between  H^^lk]  and  H3[k].  Thus  with  a  knowledge  of  one  the  other  can 
be  evaluated.  Alternatively,  since  hjn]  is  the  inverse  Fourier  transform  of  H95[k],  it  can  be 
obtained  by  an  Inverse  Discrete  Fourier  transform  (IDFT)  of  Hjilk].  Eq.  (28)  expresses  a 
relationship  between  hjn]  and  hjn],  thus  hjn]  is  known.  H^lk]  can  be  obtained  as  the  Fourier 
transform  hjn],  by  utilizing  a  discrete  Fourier  transform  (DFT)  algorithm;  This  procedure  forms 
the  basis  of  our  technique  for  the  extraction  of  a  real,  causal  time  domain  response  from  band- 
limited  complex  frequency  domain  data.  The  theoretical  development  assures  us  that  by 
computing  the  DFT’s  and  IDFT’s  the  original  real  time  sequence  will  not  lose  its  causal  nature. 


4.  INTERPOLATION/EXTRAPOLATION  OF  FREQUENCY  DOMAIN  DATA 

A  technique  to  extrapolate/interpolate  data  in  the  frequency  domain  utilizing  the  Hilbert 
Transfoirn  is  described  next.  Before  the  algorithm  is  described,  it  is  useful  to  know  something 
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about  the  available  frequency  domain  data.  Assume  that  we  have  a  complex  frequency  domain 
data  between  frequencies  f,  and  Consider  a  missing  band  between  f^  and  f,.  The  frequency 
domain  data  is  sampled  at  (n2-n,)  frequency  points  between  f^  and  f,,  and  at  (n4-n3)  points 
between  f,  and  f^.  This  is  expressed  as  a  vector 

H[n,:n,l  =  ...  0  ...  0,  H,,  ...  (38) 

It  is  now  our  objective  to  interpolate  this  missing  data  between  n2+,  and  nj.,. 

1.  The  available  bandlimited  frequency  domain  data  is  padded  with  zeros  to  ensure 

a  length  of  n  points  where  n  is  given  by  N/2  +  1,  and  N  is  [2,  4,  8,  1024, 

2048,  ...],  providing  a  sequence  of  even  length.  The  complex  data  is  now  given 
by 

H[l:n]=H[l:N/2  +  l]=[0,  0,....0,  0,...,0  ]  (39) 

2.  This  complex  sequence  is  copied  to  obtain  a  complex  consequence  of  length  N. 

This  is  done  by  flipping  the  complex  conjugate  of  the  sequence  and  appending  to  the 
given  sequence. 

H[1:N]  =  [H[l:N/2  +  l],  H*[N/2:2]]  (40) 

3.  The  complex  sequence  is  now  split  into  its  real  and  imaginary  parts. 

Hr  =  Real  [H]  (41) 

Hj  =  Iiiiflg  [H]  (42) 

4.  An  inverse  discrete  Fourier  transform  of  Hr  results  in  an  even  sequence  as  stated 
earlier. 

h^(l:N)  =  Real[IFFT(HR)]  (43) 

This  is  in  fact  the  even  part  of  the  actual  time  domain  sequence. 

5.  Before  proceeding  further,  it  is  important  to  know  that  there  are  sharp 
discontinuities  in  the  frequency  domain  signal.  In  order  to  deal  with  this  situation. 
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we  will  have  to  multiply  the  time  domain  sequence  with  a  window. 


The  Hanning  window  of  length  N  is  multiplied  with  the  time  domain  sequence. 
The  resulting  frequency  domain  sequence  will  now  be  "smeared"  or  smoothed  out  [10]. 
The  Hanning  window  is  given  by 


Hence, 


0.5  cos(27i:n)/N  0^n<N 
0  otherwise 


h^(l:N)  =  h„(l;N)  W(1:N) 


(44) 


(45) 


6.  The  odd  sequence  is  obtained  from  the  even  sequence  by  making  use  of  Eq.  (28). 

We  have 

h„(l:N)  =  [0  h^(2:N/2)  0  -  h^(N/2+2:N)]  (46) 


7.  The  discrete  Fourier  transform  of  this  odd  sequence  will  give  the  imaginary  part 
of  the  spectrum  as  stated  earlier. 

H,"®'^  =  Iinag[FFT(h^]  (47) 

8.  A  substitution  for  the  missing  points  is  made  in  the  imaginary  part  of  the  original 
sequence  using  the  sequence  obtained  in  step  (7). 


H 


Sub 


Hr'(l:nr  l)3i(ni:n2),Hi“'^(n2+l:n3-l),Hi(n3:n4),Hr>4+l:N/2  + 1)] 


(48) 


This  sequence  is  copied  to  obtain  a  sequence  of  length  N. 


H 


Subs 

I 


Hf'* [l:N/2  + 1] ,  - Hf ‘*[N/2:2] 


(49) 


which  is  an  improved  version  of  the  original  sequence  Hj. 

10.  The  inverse  discrete  Fourier  transform  of  this  sequence  will  give  us  the  odd 
sequence  again. 
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hcew 
o 


(50) 


1 1.  Again  using  Eqn.  (9),  we  get  the  modified  version  of  h^,. 

IC"  =  [he(l),  hri2:nJ2),  hJN/2  +  1),  -hnN/2+2:N)]  (51) 

12.  The  discrete  Fourier  transform  of  this  sequence  obtained  in  the  previous  step  will 
give  us  the  Real  part  of  the  spectrum  as  stated  earlier. 

=  Real[FFT(hf'')J  (52) 

13.  A  substitution  for  the  missing  points  is  made  in  the  Real  part  of  the  original 
sequence  using  the  sequence  obtained  in  step  (12). 


H„(n,:n2),Hr'(%’-l:n3-l),Hj(n,:n.^l,Hr(n4-l:N/2+l)]  (53) 


14.  This  sequence  is  copied  to  obtain  a  sequence  of  length  N. 


=  [Hr'*[1:N/2+1],  Hr“’’[N/2:2]]  (54) 

which  is  an  improved  version  of  the  original  sequence  Hr. 

15.  The  resulting  sequence  is  subject  to  an  inverse  discrete  Fourier  transform  to  obtain 
the  even  sequence. 


h„(l;N)  =  Real[lFFT(H|“’^)]  (55) 

16.  As  in  step  5,  this  time  domain  sequence  is  multiplied  with  the  Hanning  window. 

17.  Subsequent  signal  processing  are  iterations  of  steps  6  through  16. 

The  above  procedure  will  interpolate  the  missing  band  of  frequencies.  The  reconstructed 
sequence  will  now  be  the  complex  sequence  given  by 

(56) 

comparing  with  Eqn.  (38) 
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(57) 


It  is  worthwhile  to  note  that  by  making  use  of  the  Hanning  window,  although  we  have  overcome 
the  difficulties  due  to  discontinuities  at  the  ends  of  the  missing  band,  we  might  suffer  a  loss  of 
resolution.  This  is  not  a  serious  problem  and  its  effects  can  be  minimized  as  shown  in  the 
numerical  examples. 

Subsequent  sections  illustrate  this  technique. 


5.  NUMERICAL  RESULTS 

Let  us  first  consider  the  frequency  domain  data  of  a  microstrip  filter  measured  using  the 
HP  8510B  Network  Analyzer.  The  device  is  a  bandpass  filter  and  its  characteristics  are  measured 
at  415  points  from  4.31  GHz  to  7.41 5GHz.  Our  objective  now  is  to  compare  the  performance 
of  the  two  methods  (The  Cauchy  method  and  the  iterative  technique  based  on  the  Hilbert 
transform)  as  the  number  of  missing  points  are  gradually  increased.  These  missing  points  are 
created  by  deleting  portions  of  the  measured  data. 

Figure  la  shows  the  real  and  imaginary  parts  of  the  original  data.  Let  us  now  discard  40 
points  (which  is  about  10%  of  the  data)  from  200-240.  Figures  lb  shows  this  deleted  band  of 
data.  This  data  is  now  given  as  input  to  the  direct  method.  The  entire  data  set  is  not  required 
for  this  method.  Only  a  few  points  before  and  after  the  missing  band  is  stifficient.  The  program 
returns  the  interpolated  data.  Next,  this  data  is  given  as  input  to  the  iterative  method  utilizing 
the  iterative  technique  based  on  the  Hilbert  transform  presented  in  Section  4.  The  missing  points 
are  zero  padded.  Figure  Ic  compares  the  output  of  both  methods  with  the  original  real  part, 
while  Figure  Id  compares  the  corresponding  imaginary  part.  Clearly  the  reconstruction  is  quite 
accurate  using  either  techniques. 

Next,  the  number  of  deleted  points  was  increased  to  60,  i.e.,  points  200-260  were 
discarded.  The  same  procedure  was  repeated.  Figure  2a  shows  the  truncated  data.  Figure  2b 
shows  the  reconstructed  real  part  and  Figure  2c  shows  the  reconstructed  imaginary  part  of  the 
response  utilizing  both  the  techniques.  It  is  clear  from  these  figures  that  the  reconstruction 
obtained  using  the  iterative  method  based  on  Hilbert  transform  is  slightly  better  than  that  obtained 
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using  the  direct  method.  Note  that  the  amplitude  of  the  reconstructed  part  using  the  Cauchy 
method  came  out  slightly  higher  than  the  actual  amplitude. 

When  the  number  of  deleted  points  was  increased  to  80  [about  20%  of  the  data  has  been 
deleted  from  the  middle  of  the  band],  i.e.,  from  200-280,  the  iterative  method  again  proved  to 
be  better  than  the  direct  method.  But  this  time,  a  slight  modification  was  made  in  the  initial 
guess  for  the  missing  points  in  the  iterative  method.  A  straight  line  extrapolation  between  the 
ends  of  the  missing  band  was  made  in  the  initial  guess  instead  of  zero-padding  it.  Figure  3a 
shows  the  original  data  with  the  initial  guess.  Figures  3b  and  3c  are  the  reconstructed  real  and 
imaginary  parts  using  both  methods.  Clearly  the  iterative  method  gave  better  results.  But  it 
should  be  noted  that  a  better  interpolation  can  be  got  by  the  Cauchy  method  if  the  cut-off  for  the 
singular  values,  (explained  earlier  in  the  theory)  is  chosen  appropriately.  Figures  3d  and  3e  show 
the  improved  result.  The  singular  value  cut-off  was  changed  to  1  .e^'  from  the  previous  value  of 
1  .e  ’®.  However,  determining  the  cut-off  for  the  singular  values  in  practical  situations  may  not 
be  possible! 

As  the  next  example,  we  consider  the  frequency  domain  data  of  another  microstrip 
bandpass  filter  measured  at  468  points  from  4.206875  GHz  to  8.00125  GHz.  This  data  was 
incomplete.  That  is,  the  entire  data  set  was  not  considered  as  can  be  seen  in  Figure  4a,  where 
the  response  in  the  stop  band  is  not  near  zero.  Due  to  this  reason,  the  initial  guess  had  to  be 
made  as  the  straight  line  extrapolation  explained  earlier. 

Figure  4a  shows  the  real  and  imaginary  pans  of  the  original  data.  Again  we  discard  40 
points  from  200-240  as  shown  in  Figure  4b.  Since  the  missing  band  was  quite  small, 
interpolation  results  using  both  techniques  were  similar.  Figure  4c  shows  the  reconstructed  real 
parts  compared  with  the  original  real  part  and  Figure  4d  shows  the  corresponding  imaginary  part. 

Next,  60  points  were  discarded  from  200-260.  Now,  knowing  that  the  data  itself  is 
incomplete  (explained  earlier)  and  the  missing  band  is  large,  modification  in  the  initial  guess  was 
made  in  the  iterative  method.  The  initial  guess  was  again  the  straight  line  extrapolation  between 
the  ends  of  the  missing  band.  Figure  5a  shows  the  original  data  with  the  initial  guess.  Figures 
5b  and  5c  now  show  the  reconstructed  real  and  imaginary  parts  using  both  techniques.  The 
results  came  out  to  be  quite  similar.  Note  that  modifications  in  the  iterative  method  can  be  made 
just  by  observing  the  given  data  and  assessing  the  amount  of  information  provided. 
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Finally,  80  points  from  200-280  were  considered  missing.  This  meant  that  the 
interpolation  was  to  be  carried  out  for  a  very  large  missing  band  considering  that  the  data  itself 
was  incomplete.  Figure  6a  shows  the  missing  band  and  the  straight  line  initial  guess  for  the 
iterative  method. 

Figures  6b  and  6c  show  the  reconstructed  real  and  imaginary  parts  compared  with  the 
original  real  and  imaginary  parts.  Note  that  the  Cauchy  method  gave  better  results.  The 
limitation  for  the  iterative  technique  based  on  the  Flilbert  transfomi  was  that  the  given  data  itself 
was  incomplete  to  get  better  results.  So  if  the  data  set  is  incomplete  and  the  missing  band  data 
is  to  be  interpolated/extrapolated  then  the  Cauchy  method  provides  slightly  better  performance 
than  the  iterative  method  based  on  the  Hilbert  Transform. 

6.  CONCLUSION 

Both  the  direct  methods  utilizing  analytic  continuation  and  the  iterative  method  using 
properties  of  the  Hilbert  Transform  are  very  efficient  tools  for  interpolation  or  extrapolation  of 
frequency  domain  responses.  The  above  examples  demonstrate  the  implementation  of  the  two 
different  algorithms  and  compares  the  performance  of  each  method. 
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Figure  Captions 


1(a)  Original  Data  (415  points). 

1(b)  Data  with  40  missing  points. 

1(c)  Comparison  of  the  reconstruction  using  the  two  methods  (Real  Part). 

1(d)  Comparison  of  the  reconstruction  using  the  two  methods  (Imaginary  Part). 

2(a)  Data  with  60  missing  points. 

2(b)  Comparison  of  the  reconstruction  using  the  two  methods  (Real  Part). 

2(c)  Comparison  of  the  reconstruction  using  the  two  methods  (Imaginary  Part). 

3(a)  Data  with  80  missing  points  and  Non-zero  initial  guess. 

3(b)  Comparison  of  the  reconstruction  using  the  two  methods  (Real  Part).  Cutoff  in  Cauchy 
l()■'^ 

,3(c)  Comparison  of  the  reconstruction  using  the  two  methods  (Imaginary  Part).  Cutoff  in 
Cauchy  10'®. 

3(d)  Comparison  of  the  reconstruction  using  the  two  methods  (Real  Part).  Cutoff  in  Cauchy 
1()-^'. 

3(e)  Comparison  of  the  reconstruction  using  the  two  methods  (Imaginary  Part).  Cutoff  in 
Cauchy  10'^'. 

4(a)  Original  Data  (468  points). 

4(b)  Truncated  Data  with  Initial  Guess  (40  missing  points). 

4(c)  Comparison  of  the  reconstruction  using  the  two  methods  (Real  Part). 

5(a)  Truncated  Data  with  Initial  Guess  (60  missing  points). 

5(b)  Comparison  of  the  reconstruction  using  the  two  methods  (Real  Part). 

5(c)  Comparison  of  the  reconstruction  using  the  two  methods  (Imaginary  Part). 

6(a)  Truncated  Data  with  Initial  Guess  (60  missing  points). 

6(b)  Comparison  of  the  reconstruction  using  the  two  methods  (Real  Part). 

6(c)  Comparison  of  the  reconstruction  using  the  two  methods  (Imaginary  Part). 
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ORIGINAL  DATA 


Figure  1(a) 


Original  Data  (415  points) 
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RECONSTRUCTED  REAL  PART 
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RECONSTRUCTED  IMAG  PART 


Figure  1(d)  Comparision  of  the  reconstruction  using  the  two  methods  (Imaginary  Part) 
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RECONSTRUCTED  REAL  PART 


Figure  3(b)  Comparision  of  the  reconstruction  using  the  two  methods  (Real  Part)  Cutoff  in 

Cauchy  10“^® 
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Figure  3(c)  Comparision  of  the  reconstruction  using  the  two  methods  (Imaginary  Part)  Cutoff 

in  Cauchy 


25 


Figure  3(d)  Comparision  of  the  reconstruction  using  the  two  methods  (Real  Part)  Cutoff  in 

Cauchy  10“^* 
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ORIGINAL  DATA 


Figure  4(a)  Original  Data  (468  points) 
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RECONSTRUCTED  REAL  PART 
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RECONSTRUCTED  REAL  PART 


Figure  5(b)  Comparision  of  the  reconstruction  using  the  two  methods  (Real  Part) 


Figure  6(a)  Truncated  Data  with  Initial  Guess  (60  missing  points) 
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RECONSTRUCTED  REAL  PART 


Figure  6(b)  Comparision  of  the  reconstruction  using  the  two  methods  (Real  Part) 
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